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Abstract: In previous articles we outlined a subtraction scheme for regularizing doubly- 
real emission and real-virtual emission in next-to-next-to-leading order (NNLO) calcula- 
tions of jet cross sections in electron-positron annihilation. In order to find the NNLO 
correction these subtraction terms have to be integrated over the factorized unresolved 
phase space and combined with the two-loop corrections. In this paper we perform the 



integration of all one-parton unresolved subtraction terms. 
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1. Introduction 

In recent years a lot of effort has been devoted to extending the subtraction method of 
computing QCD corrections at the next-to-leading order (NLO) accuracy to the compu- 
tation of the radiative corrections at the next-to-next-to-leading order (NNLO) [1—11]. In 
particular, in Ref. [10], a subtraction scheme was defined for computing NNLO corrections 
to QCD jet cross sections to processes without coloured partons in the initial state and 
arbitrary number of massless particles (coloured or colourless) in the final state. That 
scheme can be summarized as follows. 

The NNLO correction to any m-jet cross section of processes without coloured partons 
in the initial state is a sum of three contributions, the doubly-real, the one-loop singly- 
unresolved real-virtual and the two-loop doubly-virtual terms, 

^NNLO = J dCT RR 2 j^ + J dfT HV 1 J m+1 + f da ^ J m . (L1) 

Jm+2 Jm+1 Jm 

Here the notation for the integrals indicates that the doubly-real corrections involve the 
fully-differential cross section dcr^j^ of m + 2 final-state partons, the real- virtual contri- 
bution involves the fully-differential cross section for the production of m + 1 final-state 
partons at one-loop and the doubly-virtual term is an integral of the fully-differential cross 
section for the production of m final-state partons at two-loops. The phase spaces are 
restricted by the corresponding jet functions J n that define the physical quantity. 



In d = 4 dimensions the three contributions in Eq. (|l.l| ) are separately divergent, but 
their sum is finite for infrared-safe observables. (The requirement of infrared safety implies 
certain analytic properties of the jet functions J n that are spelled out in Ref. [8].) As 
explained in Ref. [10] we first continue analytically all integrals to d = 4 — 2e dimensions 
and then rewrite Eq. ( |P| ) as 

a NNLO = / doS<'+ / rl^V°+ I d«*™>. (L.2, 



m+2 "T / {J - U m+1 T / uu m 

m+2 J m+1 



that is a sum of integrals, 

i NNLO _ f j^.RR t j RR,A 2 T f , RR,Aj T , RR,A 12 T ] \ ri Q \ 

a(J m+2 — y 1(J m+2 J ra+2 ~ ^^ m+ 2 J m ~ \&V m +2 J m+1 ~ ^^ m +2 J ™J J _ ' K 1 - 6 ) 

d«° = { [d^ +i +Jdal^} J m+1 - [do%t + ( /<+2 Al ) Al ] J m } e=o , (1-4) 



and 



m i 



a NNLO (a VV , [ \a RR,A 2 , RR,A 12 1 . [ [, RV,Ai , I f a RR.AA Ail 1 T 

d^ m u = |d< v + y [da m+2 - d<r m+ ' 2 12 j + y [do^ + { j da m+2 J J ) e= /' 

(1.5) 
each integrable in four dimensions by construction. The forms of the subtraction terms 
in Eqs. (L3) and (|1.4| ) are symbolic in the sense that each approximate cross section is 
actually a sum of many terms. The jet function depends on different momenta in each of 
those terms; the exact set of momenta for each term can be found in Refs. [10, 12]. In 
Eq. (|0|) dc7 m+ 2 1 and dcr m+ ' 2 2 are approximate cross sections that regularize the doubly- 
real emission cross section in the one- and two-parton infrared regions of the phase space, 
respectively. The double subtraction due to the overlap of these two terms is compensated 
by d<7 m+ ' 2 12 . These terms are defined in Ref. [10] explicitly, where the finiteness of d<7^+ 2 ° 
is demonstrated also numerically for the case of e + e~ — > 3 jets (m = 3). In Ref. [13], 
we computed the integral J 1 d<T m+ ' 2 1 and showed that the terms in the first bracket in 



Eq. (1.4) do not contain e poles. Nevertheless, those terms still lead to divergent integrals 



due to kinematical singularities in the one-parton unresolved parts of the phase space. In 
Ref. [12] we defined explicitly da m _^ 1 * and ( f t d<7 m+2 1 ) , that regularize the singly- 
unresolved limits of the real- virtual cross section and f. dcr, m+2 1 in turn and demonstrated 
the finiteness of the regularized cross section dcr^+ LO for the example of e + e~ — ► 3 jets. 
Thus all formulae relevant for constructing dcr^ 2 and dcr^Ki explicitly are available. 

In order to finish the definition of the subtraction scheme, one has to compute the 
integrals over the factorized one- and two-parton phase spaces, indicated in Eq. (lO). Those 
integrals have to be computed in d dimensions and the results have to be presented in the 
form of a Laurent expansion in e. According to the KLN theorem, the e poles in the 
expansions have to cancel those in the two-loop contribution da^ , leading to a finite cross 
section dcr^ NLO , that can be integrated in four dimensions. In this paper we compute the 
e-expansion of the one-particle integrals, denoted formally by f ± , that appear in Eqs. (|L~4| ) 
and (|1.5|). Once these are known, the only missing ingredient for a complete scheme 
for computing NNLO corrections is the e-expansion of the two-particle integrals, denoted 
formally by J 2 in Eq. (|1.5|). 

There are several ways to compute the one-particle integrals. If the singular integrals 
in the chosen integration variables are non-overlapping and occur in a single point in the 
integration region (which can always be mapped to the origin), then one can isolate the 
poles using standard residuum subtraction, leading to integrals of smooth functions that 
can easily be evaluated numerically. In most integrals we encounter, there are overlap- 
ping singularities in some variables and/or some variables lead to singular integrals in two 
points of the integration region. In the latter case, the integral can be written as a sum 
of two integrals with singularity in a single point, while the overlaps can be disentangled 
using sector decomposition [15], so that residuum subtraction can be applied. We have 
written a Mathematica program for the extraction of the poles employing these techniques. 
The program produces Fortran codes that can be immediately used in a Monte Carlo 



integration program such as those in the Cuba library [17]. We used the program Sec- 
torDecompisition of Ref. [16] to check our integrations. The integrated subtraction 
terms are smooth functions of a single kinematical variable (in one case three variables) 
and some parameters. 

One can also extend the method of differential equations [18], developed for com- 
puting multi-loop Feynman integrals [19], to the relevant phase-space integrations. This 
leads to e-expansions with analytic, though rather cumbersome coefficients. We have used 
this method for computing some of the singly-unresolved integrals of the approximate 
real-virtual cross section and found numerical agreement with the results of residuum sub- 
traction [20]. 

A third way to obtain the one-particle integrals is to use the Merlin-Barnes (MB) 
technique [21-23] to compute them, leading to e-expansions with analytic coefficients. The 
singly-unresolved integrals defined in this paper have also been computed via deriving MB 
representations and using the MB [24] Mathematica program for obtaining the analytic 
continuations and then the coefficients of e-expansions [25]. 

All integrals presented in this paper have been obtained by at least two independent 
computations. 

In the next section we set some basic notation used in the rest of the paper. In general 
we adopt the colour- and spin-state notation of Ref. [14] and the notation for defining the 
various cross sections introduced in Refs. [8,10,12]. In Sect. ^| we compute the integral of the 
doubly-real singly-unresolved counterterms up to 0(e 2 ) accuracy. Sects. f| and || repeat the 
same computations to 0(e°) for the real- virtual counterterms and for the counterterms of 
the integrated approximate cross section, respectively. These integrations are very similar 
to that in Sect.|3], therefore, the structure of these sections is the same as that of Sect. |3|, 
only the actual expressions differ. Sect.0 contains our conclusions. 



2. Notation 

2.1 Matrix elements 

We consider processes with coloured particles (partons) in the final states, while the initial- 
state particles are colourless (typically electron-positron annihilation into hadrons). Any 
number of additional non-coloured final-state particles is allowed, too, but they will be 
suppressed in the notation. Resolved partons in the final state are labeled by i,k,l,. . . , 
the unresolved one is denoted by r. 

We adopt the colour- and spin-state notation of Ref. [14]. In this notation the am- 
plitude for a scattering process involving the final-state momenta {p}, \-M m ({p})), is an 
abstract vector in colour and spin space, and its normalization is fixed such that the squared 



amplitude summed over colours and spins is 

\M m \ 2 = (M m \\M m ) . (2.1) 

This matrix element has the following formal loop expansion: 

\M) = \M {0) ) + \M {1) ) + . . . , (2.2) 

where \M^ ') denotes the tree-level contribution, [Ai^ 1 '} is the one-loop contribution and 
the dots stand for higher-loop contributions, which are not used in this paper. 

Colour interactions at QCD vertices are represented by associating colour charges Tj 
with the emission of a gluon from each parton i. In the colour-state notation, each vector 
\M) is a colour-singlet state, so colour conservation is simply 

(J2T 3 )\M) = 0, (2.3) 

j 

where the sum over j extends over all the external partons of the state vector \M), and 
the equation is valid order by order in the loop expansion of Eq. (2.2). 



Using the colour-state notation, we define the two-parton colour-correlated squared 
tree amplitudes as 

l^(3k)(W)r = (M^({p})\T t .T k \M^({p})) (2.4) 

and similarly the three-parton colour-correlated squared tree amplitudes, |.M^ fc «| 2 for i, 
k and I being different, 

\ M< $k,l)\ 2 =Y,^c{M^\T-T b k Tt\M^) . (2.5) 

a,b,c 

The colour-charge algebra for the product (Ti) n (T k ) n = Ti-T k is: 

T v T k = T k -Ti if i^k; T] = d. (2.6) 

Here Cj is the quadratic Casimir operator in the representation of particle % and we have 
C F = T R (iV c 2 - l)/N c = (iV c 2 - l)/(2iV c ) in the fundamental and C A = 2T R N C = N c in the 
adjoint representation, i.e. we are using the customary normalization Tr = 1/2. 

2.2 Cross sections 

In writing the final results in our notation we shall use the following cross sections: 

• the Born cross section of producing m final-state partons, 

da*=M £ d <UM)-^-|A^)(M)! 2 , (2.7) 



{m} 



s 



{m} 



the Born cross section of producing m + 1 final-state partons, called real correction, 
dal +1 =M J2 ^ m+ i(W)^— l-M^xdri)! 2 , (2.8) 

and the one-loop correction d<r^ to the m-parton Born cross section, called virtual 
correction, 

dal=M J2d4> m ({p})^—2Re(M^({p})\M^({p})), (2.9) 



{m} 



s. 



{m\ 



where N includes all QCD-independent factors and d(ft n ({p}) is the d-dimensional phase 
space for n outgoing particles with momenta {p} = {pi, . . . ,p n } and total momentum Q, 

dMPi,---,Pn;Q) = flj^p [ S + (ph(27r) d 5^ [q-J^pA . (2.10) 

The symbol Y2t n } denotes summation over the different subprocesses and S{ n y is the Bose 
symmetry factor for identical particles in the final state. 



3. Integrals of the doubly-real singly-unresolved counterterms 



The doubly-real singly- unresolved approximate cross section times the jet function reads 



d(T m+2 1 J m+1 = N J2 d( J ) rn+2({p})-^ 



E 



2 



{m+2} ^ m+2 > 

E^' 0) (M)4 + i({p} w ) 

4°'°\{p}) - Y, c ir4 m ({p})) J m+ i({p} {r) ) ■ (3-1) 

The precise meaning of the counterterms C\ r , Sr and C ir Sr ' and the definition of the 
tilded sets of momenta {p}^ tr ' and {p}^ r ' (arguments of the jet functions) were spelled out 
explicitly in Ref. [10]. The integral of this approximate cross section over the one-particle 
factorized phase space was presented in Ref. [13]. Nevertheless, we recompute it here for 
two reasons. On the one hand, we use slightly generalized subtraction terms, while on the 
other hand in a NNLO computation we need the expansion of the integrals up to 0(e 2 ) 
accuracy. In the following, we recall the definitions of the subtraction terms only to the 
extent needed to compute their integrals over the factorized phase spaces. 



3.1 Integral of the collinear counterterm 

The collinear counterterm reads 

4 0,0) (M) = ^^A« fc ^<AC 1 ({a^ ) )ip)3 P (^ P> «r,i,fc^Hc)iA^ 1 ({ia (ir) )) , 9 n , 

s ir ("J-^j 

x(l-a tr ) 2d °- 2m ^Q(a -a ir ), 

where Si T = 2pi-p r , P\ I is the Altarelli-Parisi splitting kernel in d = 4 — 2e dimensions for 
the splitting process /—►/» + /r and «o £ (0, 1]. The momentum fractions are defined as 

UiQ j J/rQ /„ „n 

and rj j = , (3.3) 



2/(ir)Q V{ir)Q 

where y iQ = s jQ /Q 2 = 2p y Q/Q 2 , (j = i, r etc.) and y (ir . )Q = y iQ + y rQ . With this 
definition, we clearly have z% r + z r ^ = 1. In order to render the integrated counterterm 
m-independent (see Eq. ( |3.6[ ) below) and to reduce the CPU time necessary for the numer- 
ical integration, we have included two harmless factors in the second line of Eq. (|3.2| ) as 
compared to the original definitions in Ref. [10]. We give a detailed discussion of these 
modifications in Appendix [A]. 



The matrix elements on the right hand side of Eq. fl3,2| ) are evaluated with the m + 1 
momenta {p}( %r > = {pi, . . . ,pi r , . . . ,p m+ 2J, that are defined by a specific mapping, 

M ^ {P} {ir) , (3.4) 

of the original m + 2 momenta {p} . This mapping leads to an exact factorization of the 
(m + 2)-particle phase space such that we have 

d</> m+2 ({p};Q) = d<f> m+1 ({p}^;Q) [dpf^ +1 ( Pr ,p ir ;Q)} . (3.5) 

We write the result of integrating the collinear subtraction term as given in Eq. ( |3.2j ) 
over the factorized phase space [dp^ +1 ] in the following form: 

/^(M) = ^S e (£j^{y~ Q ;e)Tl\M^ + .{{P}^)? , (3.6) 

where 

« __ (47T) 6 

(27r) d - 3 r(l-e)' y ' 

y-^n — s 1^qIQ 2 = ^Pir'Q/Q 2 and the arguments of the function C\ r (y~ ; e) indicate that 
this function is independent of m. In order to lighten the notation throughout the paper 
we do not show the explicit dependence of C- on ao and d,Q. 



The factorized phase space measure in Eq. fl3.5|) may be written in several equivalent 
ways. In this paper choose the form of a convolution: 

/*! 5 — 

[dPlfm+l(Pr,Pir;Q)] = / da ir (1 - a ir ) 2m ( 1_£ ) -1 -^ dfc(Pi,Pr;P(ir)) , (3.8) 




m + 2 



c„. 




(ir) ^^ r 



m + 2 



Figure 1: Graphical representation of the collinear momentum mapping and the implied phase 
space factorization. 



where pf. •. = (1 — Ui r )p^ r + a^Q 11 . The collinear momentum mapping and the implied 
factorization of the phase space measure are represented graphically in Fig. |j. The picture 
on the left shows the (m + 2)-particle phase space d(/> m +2({p}; Q), where in the circle we 
have indicated the number of momenta. The picture on the right corresponds to Eqs. ( |3.5| ) 
and (|3.8| ): the two circles represent the (m + l)-particle phase space d4> m +i({Pr ir '', Q) and 
the two-particle phase space d(f>2(Pi,Pr',P(ir)) respectively, while the symbol ® stands for 
the convolution over ai r , as precisely defined in Eq. (|3~lj|). 

With the chosen form of the factorized phase-space measure in Eq. (|3.8| ) , we can express 



the functions C ir as 



C ir (Vivo'' 6 ) 



(4tt)< 



S r 



(Q 



2\e 



(10 



dafl 



a 



,2do-l irQ 

2vr 



(0) 



(3.9) 



d4>2(pi,Pr;P(ir)) — Pfj r (Zi,r,Z r ,i; 



In writing the right hand side of Eq. (3.9) we have used that because fejLir as defined in 
Ref. [12] is orthogonal to pi r , the spin correlations generally present in Eq. (|3.2[ ) vanish 
after azimuthal integration. Therefore, we may replace the splitting kernels Pi i by their 

azimuthally averaged counterparts P\ I when integrating the subtraction term over the 
factorized phase space. These spin-averaged splitting kernels depend, in general, on Zj >r 
and z r> i, with the constraint 

Zi,r + Zr,i = 1 , (3-10) 

and are listed in Appendix |^. The phase-space measure is symmetric under the i <-> r 
interchange, thus integrals containing integrands that are linear combinations of ratios of 
the momentum fractions, as in Eq. ([3,9|), can be expressed as linear combinations of integrals 
with integrands depending only on z^. Therefore, we need to evaluate the following 
integrals: 

A- 



(47r)2 (Q 2 T f V, (1 - C^) 2 *- 1 S ~f^ fd<t> 2 (Pi,Pr;P(ir)) ^ , (3-11) 



Sr 



for the values k = — 1, 0, 1, 2. As explained in the introduction, in this paper we use sector 
decomposition and residuum subtraction to compute the e-expansion of Eq. ( 3.11|) . It is 
also possible to evaluate these integrals analytically. The details are given in Ref. [20]. 



3.2 Integrals of the soft-type counter-terms 

We refer to the soft and soft-collinear counterterms together as soft-type because they both 
use the momentum mapping appropriate to the soft subtraction term. We define these 
counterterms as follows (note again the harmless factors as compared to the definitions in 
Ref. [10]) 



x(i-yrQ) d '°- m{1 - €) Q(yo-yr Q ), (3.12) 



C ir Sr\{ P }) = SW^^Tf \M^ +1 ({p} {r) )\ 2 

$ir %r/i 
x(l-yrQ) d ^ m(1 - e) e(y -yr Q ), (3.13) 

S lk (r) = ^, (3.14) 

SirSkr 

and yo £ (0,1]- The two-parton colour-correlated squared matrix element l-Mrn+wj k \\ 2 
in Eq. ( 3.12; ) is defined in Eq. ( |2.4j ). The momentum fractions Zj jT . and z r ,i were defined 



where 



in Eq. ( |3.3j ). The matrix elements on the right hand sides of Eqs. ( |3.12j ) and ( |3 . 1 3| ) are 
evaluated with the m + 1 momenta {ftY r > = {pi, ■ ■ . ,p m +2} (Pr is missing from the set), 
obtained from the original m + 2 momenta, {p} , by a mapping, 

M -^ {P} {r) ■ (3-15) 

This mapping leads to an exact factorization of the (m + 2)-particle phase space, 

d</> m+2 ({p};Q) = d0 m+1 ({p}^;Q) [dp^ m+1 (p r ;Q)} . (3.16) 

We write the result of integrating the soft-type subtraction terms over the factorized 
phase space [dp^ +1 ] as follows: 



f4°>°H{P}) = %s t (iL) EEs!?(%^)i<;(„„(® (r) )i 2 . (3-i7) 

lc ir Sr\{ P }) = ^S e {£)c^(e)T 2 \M%{{p) ir) )\ 2 ■ (3-18) 

The functions S^ fc (Y^, Q ;e) and CS^ '(e) are again independent of m as the arguments 
indicate and as before, we lighten the notation throughout by not showing the explicit 
dependence of the integrated subtraction terms on yo and d' . Also, S 4 - fc (Y^ Q ;e) only 
depends on the combination of variables 

%Q = 7^t- ■ (3-19) 

UiQUkQ 



As in the collinear case, the one-particle factorized phase space may also be written in 
the form of a convolution 



/•l 

[dPlL + l(Pr;Q)] = / dy rQ {l-y rQ 
Jo 



m(i-f)-lV_ ( 

2vr ( 



i(Pr,K;Q) 



(3.20) 



where the timelike momentum K is massive with K 2 = (1 — y T Q)Q 2 ■ We show the soft 
momentum mapping and the implied phase space factorization in Fig. 0. The picture on 




m + 2 



S,. 




m + 2 



Figure 2: Graphical representation of the soft momentum mapping and the implied phase space 
factorization. 



the left shows again the (m+2)-particle phase space dcj) m +2({p}', Q), while the picture on the 
right corresponds to Eqs. ( p. 16 ) and ( 3.20[ ): the two circles represent the two-particle phase 
space dfaipr , K ', Q) and the (m + l)-particle phase space d(f) m +i({p}( r > ; Q) respectively. 
The symbol (8> stands for the convolution over y r Q as defined in Eq. (|3,20|). 



Using the definition of the subtraction terms in Eqs. (3.12) and ( |3.13j ) and the chosen 
form of the factorized phase space, we find that the functions S^. (l^. Q ;e) and CS^ ^(e) 
can be written as 



q(0W . A _ ( 47r ) 



CS( 0) (e 



(4^) 2 



dy(l-y) d o- ll ±-dMPr,K;Q)-S tk (r) 



2\e 



m 



2vr 



( \ !){ l-yfo-^dMPr,K;Q)-^ 

S{ r %r.i 



(3.21) 
(3.22) 



The computation of these integrals is fairly straightforward using energy and angle vari- 
ables. The details can be found in Ref. [20]. 



3.3 The integrated approximate cross section 



DD A 

We are now in a position to compute the integral of d<7 m+ ' 2 x as given in Eq. ( |3.1| ) over 
the one-particle factorized phase space. In order to evaluate J ± do" m+ ' 2 1 we first use the 



phase-space factorization properties of Eqs. ( |3.5[ ) and (3.16), then perform the integration 



10 



to obtain 






1 a s 



{m+2} 



S{m+2} 27r VQ 



^e ( 7^ 



2 \ e 






+ E s ^ 0) (%,Q; e )l^i;M)(W)l 2 ^+i(W) 



CSM(e)T*\M { Zi({p})\ 2 Jm + i({p}) 



(3.23) 



This result is not in the form of an (m + l)-parton configuration times a factor. In order 
to rewrite Eq. ( |3.23| ) in such a form, we need to perform the counting of symmetry factors 
for going from m + 1 tom + 2 partons. This was done in Appendix B of Ref. [13] for going 
from m to m + 1 partons and the results can readily be taken over by setting m — ► m + 1. 

The final result for the integral of the doubly-real singly- unresolved approximate cross 
section can be written as [13] 



d^+' 2 Al =d^ +1 ®/ (0) (W m+ i;e), 



(3.24) 



where the insertion operator reads 



I { °H{p}r 



l+ i;ej 



o 



2ttMq 2 



J]) C|°W e) T? + J] S^(Y lktQ ; e) T v T k 

kj^i 



In Eq. ( 3.25p we introduced the functions 



Ci 0) = C<°> , C<°> = leg + nfCi? , §£> = si? + csw 



w 



'W 



(3.25) 



(3.26) 



with C\% sjf and CS (0) defined in Eqns. Q, (g33J) and ( gjgg ), respectively* In this 
paper rtf denotes the number of light flavours, which we set to rtf = 5 in all numerical 
results. Upon expanding these functions in e we find that the coefficients of the poles in 
the expansions are independent of the cut parameters «o an d yo as wen as the exponents 
do and d' and they read 



C<°>(*;e) 



I + I(f _21nx)+0(e°^ 



C<°>( 



e \2 
1 /ll 



x;e 



2 T R 

+ 7 It " T^ci 



21nxJ + 0(e°) 



Sj° } (^;e) = -lnl- + 0(e° 



(3.27) 
(3.28) 
(3.29) 



*Note a slight rearrangement of terms with respect to Ref. [13]: here we find it more convenient to use 
colour-conservation (Eq. (2.3)) to combine CS 1 - ' with S\ k into S>\ k . Then of course C^ is defined without 
including CS (0) . 



11 



For the coefficients with non-negative powers of e we obtain integral representations that 
can be evaluated numerically. The results for four values of ao or hq = 1, 0.3, 0.1 and 0.03 
with fixed values of do = d' = 3 — 3e are presented in Figs.||| and |||. The dependence on 
ao in the collinear functions is hardly visible. Note also that the small- x behaviour of the 
collinear functions is dominated by the logarithmic terms that are the same in the quark 
and gluon functions, therefore, the expansion coefficients of C q (x;e) and C g (x;e) are 
very similar. 



Using colour-conservation (Eq. ( |2.3D ), the definition of Y^q ( 3.19| ) and Eqs. ( ^27 )- 
(3.29) above, it is straightforward to verify that our insertion operator differs from that 
defined in Eq. (7.26) of Ref. [14] only in finite terms, i.e. 



/«*(«,,.)-£«,(£ 






+ 0(e' 



0\ 



with the usual flavour constants 



Po 



lq = 2 Cf and 7s O 



(3.30) 



(3.31) 



It follows that j x d<T m+ ' 2 1 correctly cancels all e-poles of the real- virtual cross section da™ + 



+i- 



4. Integrals of the real-virtual counterterms 

The approximate cross section that regularizes the singly-unresolved limits of the real- 
virtual cross section times the jet function is 



da m ^ 1 1 J m =Af Y^ d( Prn+i({p})-y — - 

{m+l} 

*{e[ ' 

I r L i^r 



, , , ~>{m+l} 



i^r 



+ E 



E^'iiiiiip) 



(ir)^ 



i=£r 






(4.1) 



The subtraction terms C\ r ' , Sr and C ir Sr ' (with (k, I) = (0, 1) or (1,0)) that appear 
in Eq. ([O]) above were explicitly given in Ref. [12]. Here we compute the integral of the 
approximate cross section over the factorized one-particle phase space. 
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4.1 Integrals of the collinear counterterms 

The collinear counterterms C\ r ' ({p}) and C\ r ' ({p}) are defined as 



x (1 - a ir ) M »- 2 ( m - 1 )( 1 - e )0(a o - «ir) , 



(4.2) 
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4 1,0) (M) = (8™ s ^) 2 ^crco S (7T6) 

x(MW({p}^)|Pg r (^,^,^ ir ;e)|^)({p}^)) 

x (1 - a ir ) M °- 2(m - 1)(1 - e) e(a - a ir ) . (4.3) 

The Pf.fX z i,r,z rt i,k± > i jr ;e) kernels are the same tree-level Altarelli-Parisi splitting func- 
tions that appear in Eq. ( |3.2| ), while the P* I (<Zi, r , z r ^, k±^ :r ; e) kernels are the one- loop gen- 
eralizations of the Altarelli-Parisi splitting functions as spelled out explicitly in Ref. [12]. 
The m momenta, {p}^ ir \ entering the matrix elements on the right hand sides of Eqs. ( |4,2j ) 
and ( |4.3| ) are obtained from the original m + 1 momenta, {p} , by the specific momentum 
mapping of Eq. ( |3.4| ) 

M ^ {p} (ir) ■ (4.4) 

This mapping leads to an exact factorization of the (m + l)-particle phase space. The form 



of this factorization is identical to the one discussed in Sect. 3J. In particular, Eqs. ( |3.5| ) 



and ( |3.8D remain valid after making the replacement m — ► m — 1. 

We write the integrals of the collinear subtraction terms over the factorized phase space 



as 

,2 \ e 



I 



[CS' 1} ({P}) = ^S e (iL) ct\y~ Q ;e)Ti2Re(M^m^)\\Mg\{p}^)) (4.5) 



and 



with 



Cr co S (7re)c£ ) (y^ Q ; e )T? r |^)({p}(-))| 2 , (4.6) 



cr= I m + e^l-e) 
r (4vr) 2 -^ r(l-2e) ' K ' 



Note that 

\2 



^ ««■« - m$t^ ^-v- **■ + & + ^ • «-» 



In writing Eqs. ( |4.5| ) and ( [4.6[ ) we used that by the usual arguments, the spin correlations 
generally present in Eqs. Q4.2| ) and ( ^.3| ) vanish after azimuthal integration so that we can 
replace the splitting functions by their azimuthally averaged counterparts when computing 



,(0) 



the integrals. The C ir (y~ ;e) function appearing in Eq. ( |4.5| ) was computed in Sect. 3T 
while C\l>(y~; ; e) in Eq. (O) is given by 



J ir \y irQi*-) g 



fOtQ ST" 1 1 

x / da(l-a) 2d °- 1 ^^dcj)2(pi,Pr]P(ir))^^PfUzi,r,z r ,i]e)— T . 
Jo 2vr s - r +e JlIr Tf r 



(4.9) 
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We note that C ir (y-r^ ;e) is independent of m as the arguments of the function indicate 
and its explicit oq and do dependence is suppressed in the notation as usual. 

Inspecting the actual form of the one-loop Altarem-Parisi splitting functions as recalled 
in Appendix [B| and using the symmetry property of the factorized phase space under the 
interchange i <-> r, we find that C\ r (y~-e) can be expressed as a linear combination of 
the integrals 



^(Q 2 )- jTWd-a*) 2 *- 1 ^ 



d<f> 2 (p 



^k+8t 



,(±)/ 



hPr;P(ir))^+T9i {Zr,i), (4.10) 



for k = —1,0,1,2 and the values of 5 and g\ as given in Table H. Here we compute 



5 


Function 


«^M 





9A 


1 


Tl 


s f> 


(1 - z) ±e 





& 


(l-z) ± * 2 F 1 (±e,±e,l±e,z) 


±1 


«g° 


2 F 1 (±e,±e,l±e,l-z) 



Table 1: The values of S and g\ (z r ) at which Eq. ( 4.10 ) needs to be evaluated 



the e-expansion of these integrals using sector decomposition and residuum subtraction. 
Analytic results are presented in Refs. [20,25]. 



4.2 Integrals of the soft-type counter-terms 



The soft and soft-collinear counterterms appearing in Eq. ( fl.l| ) are 



4°' 1} ({P}) = -8vra sM 2 ^^^ fc (r)2Re(^)({p}W)|T r T fe |^«({p}M)} 



(4.11) 



C. Si ' 1 ) = 87ra s u 2e 



1 2js 



Sj r %r 



Tl2Re(M$({p}V)\\M£\{p}V)) 



i,r rr2 



x(i-yrQ) d °~ {m - 1){1 - e) e(m-yr Q ), 



(4.12) 
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4 1,0) ({P» = (8vra s ^) 2 c r cos(Tre) £ £ ~«S ife (r) 



i fc^i 



Ca— -^ — - I t^O") + 



+ 2- 



e 2 sin(-7re) \2 
7T 1 



e cos 7re 



2e (47r) 2 cr ^ 2e cos(7re) 
E(^(r)) £ |< (w) ({p} W )| 2 } 



,(0) 



SlW^I 2 



\M^ m (W r) )\ 



l^i,k 



x(l-y rQ ) d °- im ~ 1){1 ~ e) @(vo-y rQ ) 



(4.13) 



CtrSpM 



37ra s yU ) crcos(7re) — T„- 



67V Zr> 



C A 



1 7re / 1 z. 



_, ] 0o S t 1 

e 2 sin(7re) \si r z rj i ) 2e (47r) 2 cr /i 2e cos(-7re) 



IA^WpW 



x(l-yr Q ) d »- (m ~ 1)(1 - £) e(y -yrQ), 



(4.14) 



with 



A) = y C A - -T R n f - -T R n s . 



(4.15) 



In QCD, for the number of scalars we have n s = 0. The three-parton colour-correlated 



squared matrix element l-M^-fi fc nl 2 appearing in Eq. ( 4,13| ) is defined in Eq. fl2.5| ). The 



m momenta {pp r ' that enter the squared matrix elements on the right hand sides of 
Eqs. fl4,ll|) - fl4.14j) are obtained by applying the momentum mapping of Eq. ( |3.15|) to the 
original set of m + 1 momenta {p} , 



M -^ {P} ir) 



(4.16) 



As discussed in Sect. ^2, this mapping leads to an exact factorization of the (m+l)-particle 
phase space. Indeed, Eqs. ( |3.16 ) and ( |3.20| ) remain valid after making the replacement 
m—nn—1. 



Integrating the soft-type subtraction terms of Eqs. ( [4.1l| )-( |4.14| ) over the factorized 
phase space, we can write the results as 



lsr\{p}) 



(1 \ e 
K 
Q 2 
3 (o) 



(4.17) 



EEC^i^^^^lTrr^^f))), 



i k^i 



j^C ir S^\{p}) = ^S t (jpj CS^(e)T 2 2Re(M^({p}^)\\M^({p}^)), (4.18) 
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^• 0) (M) 






^' z (4vr) 2 



S e 



cr cos(7re) 



f(°). JSZ\(r)\\2 



X EE sS } (% jQ ; e )I^S i>fe) (-b~} . 

+ £ s !2(%<3.%.%Q^)i<w)(© (r) )i 2 

l^i,k 



(4.19) 



C fr 5^' )({p}) 



(2 \ e 
Q 2 



(4tt) 



-cr cos 



(^)CS«( e )T 2 |A^({p}M)| 2 . ( 4. 2 o) 



The notation anticipates that S ik , S ikl and CS^ 1 ** are m-independent. (As usual the ex- 
plicit yo and d' dependences are not indicated.) The tree-level functions S ik (^q;c) and 
CS^ '(e) are given in Eqs. ( 3.21 ) and (|3.22| ) respectively, while their one-loop counterparts 
read 

^ik ( Y ik,Q' e ) = 



r 1 7re (47r)2 tcP 

e z sm(7re) b e j 



</l) 



(Q 2 ) 2e I' dy(l-yf^ l ^dMPr,K;Q) ( -S llA r 



2vr 



1+e 



' 2e (4vr) 2 c r 



,, j cos(7re) 



-i 



s(o) 



^ik 0%,Q) e ) i 



2 
CS (1) (e) 



(4.21) 
(4.22) 



^ 1 (^) 2 ( ^2e ''"" 
e cos(7re) 5 6 ./,, 



(Q-) 2 ' / fl.'y(.i-.v)" ; '- 1 ^d^ 2 (p r ,if;Q)i5 ifc (r)(^, / (/-)) 



\<ft- 



I ™ (4^) 2 ^ 2n2£ /•"" 

/(i 



' e 2 sin(7re) S* e 



-g A ^ ■ 7 77' (Q*) 2e / dy(l-y) d o-^d^ 2 ( Pr .,^;Q)2(— ^ 



_iQ_ 2 

2vr' 
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1+e 



go g e 

2e (4vr) 2 c r 



, , cos(7re) 



CS (0) (e) 



(4.23) 



These integrals can also be evaluated using energy and angle variables The explicit com- 
putation is performed in Ref. [20]. 



4.3 The integrated approximate cross section 

The computation of L da m ^ 1 (including the counting of symmetry factors) proceeds 
along the same lines as that in Sect. |3.3| . The final result for the integral of the real-virtual 
singly-unresolved approximate cross section can be written as (cf. with Eq. fl3.24|) ) 



/ 



d^f 1 = dal ® l(°)(M m! e) + dal ® I (1) ({ P } m , e) , 



(4.24) 
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where the insertion operator J' ' is given in Eqs. fl3.25|) -( |3,26D , while /' ' reads 



I {1 \{p}r 



2ir 



S r 



Q 2 



r 2 (i + 6)r 4 (i- e ) 

r(l + 2e)r 2 (l-2e) 

(1)/ 



E C * (1 W «0 T\ + E S«(Y^ Q ; e) T % -T k 

+ LL S ikl ( Y ik,Q> Y il,Q> Y kl,Q' £ ) 2_, fabcTiT k T{ 
kj^i ly^i,k a,b,c 



(4.25) 



In Eq. ( 4.25J ) we introduced the functions 

c (i) = c (i) 

9 99 ' 



C« = ic«+n f cS, f = S W + CS« 



yq 



(4.26) 



with C^ } , SJJ and CS (1) defined in Eqs. flO|), ( ggg ) and (|4T23|) , respectively. The terms 
proportional to (3q in the one- loop functions (superscript (1)) are proportional to the corre- 
sponding tree- level contributions (same function with superscript (0)), with proportionality 
factor 

A) 



2e 



Q2 ) cos(;T ° 



(4.27) 



and represent the effect of one-loop UV renormalization. In the following we present results 
for the unrenormalized functions CW m , C^^ and S^. bare , obtained by setting /3q = 0. 

Expanding these functions in e, we find that the coefficient of the two leading poles 
are independent of the cut parameters ao and yo as we U as the exponents do and d' : 



C g,bare( X ; £ ) 



C 



(1) 
g, bare 



x;e) 



1 

1 



C A 



1 



e c 



lnx C A + 0(e^ 2 ) 



SjfcharcV* ! £ ) 



S^(n,Y 2 ,Y 3 ;e)=^ + 0(e-- 



r l f 11 i _l X Tr 
"4? CA -?(T2- lnX+ 3 nf ^ 

-^c A ^ + o( e - 2 ), 

3vr 



2 TrCf 

3 nf ^r 



(4.28) 

C A + 0(e- 2 ), (4.29) 

(4.30) 
(4.31) 



For the remaining coefficients we obtain integral representations which can be evaluated 
numerically. The results for four values of ao or yo = 1) 0.3, 0.1 and 0.03 with fixed 
d' = 3 — 3e are presented in Figs.|5|-p|. As in the case of C- functions, the 



values of do — u, 

dependence on ao in the collinear functions is hardly visible and the expansion coefficients of 
C^ l arc {x; s) and C bare (x; e) are very similar. The function S ikl depends on three variables 
Y\, Yi and I3. In the plots we only show representative results obtained by fixing two 
variables at 0.1 and varying the third. 
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5. Integrals of the counterterms for the integrated approximate cross sec- 
tion 



The singly-unresolved approximate cross section regularizing the integrated approximate 
cross section J t dcr, m+ ' 2 * times the jet function reads 



(/d^+' 2 Al ) Al J m =A^ £ d^ m+1 ({p})-^— 
KJl ' {^1} ^ m+1 > 



E 



E^r'(M)4(M (ir) 



i=£r 



+ 4 W) ({P» " T, C ir4 W H{p}) Jm({p} (r) ) 



rj^i 



+ £ 



E;C (0 ' 0> (W)^({?} 



(ir)^ 



i^r 



+ s?* m m) - £c^ x(0 ' 0) (M) J m ({p} w ) 



(5.1) 



The precise meaning of each subtraction term appearing above was spelled out explicitly 
in Ref. [12]. 

5.1 Integrals of the collinear counterterms 



The collinear counterterms are 



7(0,0®/), 



1 



Sir 



and 



x (1 - airf^'^-^-^Biao - a ir ) (5.2) 



, fix (0,0), 
'ir 



C2* y " u> ({p}) = 87ra s /i 2e — RiriVir, Zi,rV~: n , z r ,iy~ n ;e) 



irQ' 1 r ' 2y ir( 

x{l-a ir ) 2d °- 2 ^~ l ^ l -^Q(a -a ir ), 



(5.3) 



where P\ I are again the tree-level Altarelli-Parisi kernels and the momentum mapping 
used to define the momenta entering the matrix elements on the right hand sides and 



the corresponding phase space factorization are the same as in Sect. 4.1. The function 
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^HryViri z i,rV j r Q> ^ r >*2/jrQ' ^' r ^adS 



RiriVir, Zi, r y ir Q, z r,iy ir Q~, e) — 2 7T 1 Q 2 



cf^,^; 6) T? + C^\z r ,y~ Q , e) T 2 - C^ r) (y~ Q , e) T% 



(5.4) 



+ (T? r -T?-Tj!)S 



2\q(0) / 2/ir 



^i.r^r.iV' 



,e 



Note that in order to lighten the notation, we do not explicitly indicate the dependence of 
Rj r on ao, do, yo and d' . 

As usual, we write the integrals of the collinear subtraction terms over the factorized 
phase space in the following form 



r (0,0®J) _ «s q ( ^X C {Q) (n~ -AT 2 



x(M^H{p} {ir) )\I m ({p} {ir) -,e)\M^({p}^)) 



(5.5) 



and 



C 



Rx(0,0) 






(4tt) 



; cf (0) fe o ;e)Tf r |^)({p}W)| 2 . (5.6) 



The functions C£ J and J (0) in Eq. Q are defined in Eq. Q and Eqs. (^lEj)-(¥M), 



while C ir reads 



'in 



cfr Xi °\y^o^) = (Q 2 y / dall-aj^-^d^fe,^;^)) 



2tt 



(0) 



x i>-jr(,yir) z i,r]J ^ r Qi z r,iU i r n'i e J fif r \ Z i^^ z r,i'i e J ^,2 



(5.7) 



where 



R; 



2~7T 



^Sj^r 



Q 2 



R/ 



(5.8) 



As the rest of the integrated collinear functions, C ir too is independent of m as the 
notation suggests. Also in keeping with our conventions, we do not indicate the explicit 
dependence on the cut parameters «o and yo or the exponents do and d' . The analytic 
evaluation of the integrals is performed in Ref . [25] . 
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5.2 Integrals of the soft-type counter-terms 



The soft-type counterterms read 

i k=£i 

x (M^H{p}^)\{I({p}^;e),T v T k }\M^({p}^)) 

x(l-y rQ fo-^-^-^e(y -y rQ ), 

C^r^ttp}) = S™ s ^- 2 -^n (^)(MM)|/({p}W;e)|MM(MW)) 
x(i-yr Q ) d «- (m - 1)(1 - e) e(y -yrQ), 



(5.9) 



(5.10) 



and 



5 i?x(o,o) (W) = -8na s n 2e Y,Y,l S ^ r ^(y 



iki Vir,ykr, ViQ, UkQ, UrQ', e) 



i k^i 



x \M®. (i k)({p} {r) )\ 2 ( l ~ yrQ) d '°- {m - 1){1 - €) e(y - y rQ ) 



m;(i,k) ' 



(5.11) 



C ir S? x ^°\{p})=8ira sf i 



2c 



1 2 ^,r T 2 c ^g (\^_ 
Sir Zf A Z7T y i^/ 



C(°)( yrQ ;6)-S; 



(0) / Vir 



ViQVrC 



X |-M^({p} (r) )| 2 (l _ y rQ fo-(m- D(l-Oe (yo _ y r( 



(5.12) 



The momenta which enter the matrix elements on the right hand sides and the corre- 
sponding phase space factorization are the same as in Sect. 4.2. In Eq. ( |5.11 ) the function 

Rik,r(llik, yiriVkr, ViQ, VkQiVrQ] e) IS 

,2 \ e 



Rik,r{yik,yir,llkr,yiQ,ykQ,yrQ;e) = C*A — S t { , 



•27T^VQ2 



Cf(y rQ ;e) + sS 



cj"" I ■■'" I o(0) ( Vkr . 



(o) ( Vik. 



(0) / _J/ir_ . 



(5.13) 



.yiQykQ j \yiQy r Q 

where again we do not indicate the dependence of Rifc,r on ao, do, yo and d' explicitly. 

After integrating the soft-type subtraction terms over the unresolved phase space we 
can write the results as 

,2 \ e 



x ^(^({p} (r) )IU({p} (r) ;6),T r T fe }|^)({p}M)) , (5.14) 

J^Sr^ttp}) = ^S e (|^ CS<°>(e)T? 

x (M^({p}^)\I({p} {r) ;e)\M^({p}^)}, (5.15) 
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s^°\{ P }) 



rf^lM) 



27rMQ2 



'^' z (4tt) 2 



(0) 



SV 



EEs2 x(0) (% Q ; e )l< (i , fc) (W 



lW^i 2 



i fc^i 



(5.16) 



"s, //x 2 Vl 2 (47r) 2 
27r^VQ 2 /J ^ 



CS fl x(°)(y^ o;e )T 2 |^)({p}M)| 2 . (5.17) 



The S ik and CS^ ' functions are defined in Eqs. ( 3.21 ) and (|3.22| ) respectively, while for 



S-, and CS Rx (°> we need to compute the integrals 



-,i?x(0). 



Wo 



,-iQ 2 



i 



Sr W (% Q ;e) = -(QTy o dy(l-yr«-^dMPr,K;Q)-S tk (r) 

X Ca Rik,r{]Jik, UiriVkr, UiQiVkQ, VrQ] e) 



(5.18) 



and 



CS Rx (°) (e) = (Q 2 r jf dy (1 - yfo-^dMPr, K; Q) 



1 2z, 



Sir %r.i 



xC A 



ViQVrQ 



In Eq. (|5~lp 



^Hk,', 



Ca 2^ VQ 1 



2 \ e 



-1 



Rife,: 



(5.19) 



(5.20) 



The soft functions S ifc and CS ^ are again independent of m as the notation suggests 
but do depend on the cut parameters ao and yo as well as the exponents do and d' . The 
dependence on the later four parameters is suppressed in the notation as usual. The 
analytic evaluation of the integrals is performed in Ref . [25] . 



5.3 The integrated approximate cross section 

The computation of f ± ( f ± do" m+ ' 2 1 J (including the counting of symmetry factors) pro- 
ceeds along the same lines as that in Sect. |3.3| . The final result for the integral of the 
iterated singly-unresolved approximate cross section can be written as (cf. with Eq. ( |3.24| )) 



l x ji 



da^) M =^ 



l{l^({p} m ,e),I^({p} m ,e)}+I R ^({p} m ,e) 



, (5-21) 



where the insertion operator 1^ ' is given in Eqs. ( |3.25| )-( |3.26D , while i" *' ' reads 



r flx(0) 



({p} m ;e) 



(o \ en 2 
Q 2 



£ cf (0) (^; e )T 2 + ^sf( )( %Q ; e )T r T fe 

k^i 

(5.22 
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In Eq. ( 5.221 ") we introduced the functions 



p.Rx(0) 



^iJx(O) 



(47T) 2 



(4vr) 2 



pflx(O) 



pcRx(O) 
9 



(4tt) 5 



S r 



I C i?x(0) , c «x(0)\ 
2 99 ^ f 99 I ' 



,flx(0) + QgRx(0)\ 



(5.23) 



In this paper we evaluate the e-expansion of these functions using sector decomposition. 
We are able to find the coefficients of the two leading poles analytically (the coefficient of 
the 1/e 4 pole is zero in each case): 



CfC°)(z; £ ) = i 



- - In a? - -n f — + £(y , d' ) 



C flx(0) ( 



x;e) 



11 , 2 T R C F 

y-lnx- -nf—— + Y,(y ,d Q 



C A + 0(e- 2 ) , 
C A + 0(e- 2 ) 



and 



s2 x(0) (y ;£ ) = lc A i ny — 



+ 0(e" 



(5.24) 
(5.25) 

(5.26) 



In Eqs. ( |5,24| ) and ( 5.25| ), the function S(yo,dg) depends on the cut parameter yo and on 
the exponent d' Q . If we set d' Q = D' + d^e, where D' is an integer (greater than 2, see 



Appendix A. 3 ) and d[ is real, then we find 



D> 



£( y oX)=in yo -x:-~ (1 ~ yo)fc 



fc=i 



A- 



(5.27) 



For the remaining coefficients we obtain integral representations which we can evaluate 
numerically. The results for oq = Vo = 1> 0.3, 0.1 and 0.03 with fixed values of do = d' = 



3 — 3e are presented in Figs. 10 and O. 



6. Conclusions 

In this paper we computed the integrals over the phase space of the unresolved parton of 
the singly-unresolved counterterms of the subtraction scheme defined for computing QCD 
jet cross sections at the NNLO accuracy in Refs. [10, 12]. The results are given in the 
form of Laurent-expansions in the regularization parameter e keeping the first five terms in 
the expansion. We computed the coefficients of the two leading poles analytically, and the 
remaining three coefficients numerically. The final forms of the integrals are combined into 
insertion operators times various cross sections with the same number of external legs as 
in the doubly-virtual cross section so their combination into a single numerical integration 
is straightforward. 

The necessary integrations were carried out using iterated sector decomposition and 
residuum subtraction. In order to check the rather cumbersome computations the same 
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Figure 10: Expansion coeffiecients of the functions C i x (x; e) with do = d' = 3 — 3e and m = 5. 
Upper row: i — q, lower row: i = g. 
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Figure 11: Expansion coeffiecients of the function Sj fc bare (Y; e) with d' Q = d = 3 — 3e and ri{ = 5 



integrals were computed by analytical techniques as well. The details of those works are 
presented in separate articles [20,25]. 

We presented the expansion coefficients of the integrated subtraction terms in the form 
of plots. These plots only serve to demonstrate that in spite of the large complexity of the 
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integrands, the resulting functions are smooth, therefore, high-precision approximations can 
be easily obtained. In an actual computation of a physical observable at NNLO accuracy 
one needs only the 0(e°) piece in its expansion, for which one can use any approximation 
that describes the function with better than the expected relative accuracy of the observable 
in the kinematical region of interest. The coefficients of the lower orders are needed only to 
check the cancellation of the epsilon poles, which has to be done only once independently of 
the observable. Since the numerical integrations that compute the integrated counterterms 
converge quickly, any resonable precision of the cancellation can be achieved. 
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A. Modified subtraction terms 

In this Appendix, we describe a simple modification of the NNLO subtraction scheme 
presented in Refs. [10,12]. The purpose of the modification is twofold: first, we wish to 
make the integrated subtraction terms independent of the number of hard partons, to. 
Secondly, to save CPU time and have a better control on the numerical calculation, we 
restrict the phase space on which the various counterterms are subtracted. This has been 
found useful in the context of NLO calculations [26-28]. 

We recall the details of our NNLO subtraction scheme only to the extent we need to 
define the modification of the singly-unresolved counterterms. For further details, we refer 
the reader to the original papers. 

A.l Modification of singly-unresolved subtraction terms 

We first consider the singly-unresolved approximate cross section d<r m+ ' 2 1 appearing in 
Eq. fll.3;). We write this term in the following symbolic form: 

d<W* 2 Al = ^ m+1 [A Pl ]Ai\M { Z2\ 2 , (A.l) 

where the singly-unresolved approximation ^il-M^^I 2 * s a sum °f collinear, soft and soft- 
collinear terms (see Eq. ( |3.1| )). The precise definition of these three terms involves the 
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specification of two momentum mappings (see Eqs. (|3.4j ) and fl3.15|) ) 

M ^ {P} {tr) , {p} -^ {P} {r) , (A.2) 

both of which lead to an exact factorization of phase space in the form 

d(j) m+2 {{p}) = d(f> m +i({p})[dpi tm+ i] . (A.3) 

The exact form of the one-particle factorized phase space [dpi i?n+ i] is irrelevant for the 
present argument, its only feature which is important for us now is that it depends on the 
number of hard partons m through factors of (1 — aj r ) 2m ' 1_e ' 1 and (1 — yrQ) m ^ 1 ~~' L ' 1 for 
the collinear and soft mappings respectively, i.e. 

[dpg +1 ] oc (1 - a %r f^-^ , (A.4) 

[drf!*] « (1 - VtqY* 1 -*- 1 (A.5) 



(see Eqs. (3.8) and ( 3.20J )). The subtraction terms, as originally defined in Ref. [10], are 



m- 



independent, therefore, the m-dependence of [dpi im +i] is carried over to an m-dependence 
of the integrated subtraction terms. However, this dependence on m enters in a rather 
cumbersome way in the integrated counterterms (see e.g. Eqs. (A. 9) and (A. 10) of Ref. [13]). 

It is then worthwhile to reshuffle the m-dependence of the integrated counterterms into 
the subtraction terms themselves, where it appears in a very straightforward and harmless 
way through factors of (1 — «j r ) and (1 — y r Q) raised to m-dependent powers multiplying 
the original subtraction terms as in Eqs. (|3.2| ), ( |3.12 ) and ( 3.13 ). These factors do not 



influence the behavior of the subtraction terms in the infrared limits because 

C ir (l - a ir ) = 1 , and S r (l - y rQ ) = 1 , (A. 6) 

where Cj r and S r are the symbolic operators to take the collinear limit of momenta p^ 
with pr, or the soft limit of momentum pp, as defined precisely in Ref. [8]. Furthermore, 
the modified soft-collinear term still cancels the overlap of the soft and collinear terms 
correctly due to 

S r (l-a ir ) = l. (A.7) 

Thus, we are free to include the factors of (1 — c%.) and (1 — y r Q) as in Eqs. fl3,2| ), Q3.12J ) 
and ( p3| ). 



The 0-functions that also appear in Eqs. ( |3.2[ ), ( |3.12| ) and ( 3.13J ) control the region of 



the (m + 2)-particle phase space over which the subtraction is non-zero such that ao = 1 
and/or yo = 1 corresponds to subtracting over the full phase space. In a NNLO computa- 
tion the (very) large number of subtraction terms and their complicated analytic structure 
makes the evaluation of the approximate cross sections rather time consuming. Constrain- 
ing the phase space over which the subtraction is non-zero can result in large gains in CPU 
time. The introduction of the cutoffs also provides a strong consistency check on the whole 
of the calculation: the final results should be independent of «o and Vo- Finally, suitably 
chosen values of the cut parameters can actually improve the numerical behavior of the 
code. 
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A. 2 Modified real-virtual subtraction terras 



The real- virtual subtraction terms da m ^ 1 and f j x do" m+ ' 2 1 J which appear in Eq. (1.4) 
read symbolically 

daZlf 1 = dcf> m [d Pl ] Ai2Re(M { Zi\\M% +1 ) , (A.8) 

and 

(l^T+^Y 1 = dcf> m [d Pl }A 1 (\M^ + x\ 2 ^ (0) (W m+1 ;e)) • (A.9) 

The singly-unresolved approximations 

^Re^+illA^i-l) and A^M^U^I^ttphe)) (A.10) 

are sums of collinear, soft and soft-collinear terms (see Eqs. (|4.1|) and (|5.1|)). The subtrac- 



tion terms in Eqs. (4.1) and ( |5.1| ) are all defined using the singly-unresolved momentum 
mappings of Eq. (|A.2|) , which now map the original set of m + 1 momenta appearing in the 
one-loop squared matrix element into a set of m momenta. The appropriate phase space 
factorization reads 

d(j> m +i({p}) = d<p m ({p})[dp ltm ] , (A.ll) 



which is of course just Eq. (A. 3) with the substitution m — ► m — 1. Again the factorized 
one-particle phase spaces carry an m-dependence through factors of (1 — c%.) and (1 — y r Q) 
raised to m-dependent powers. For the collinear and soft mappings we have respectively 

[dpgj] oc (1 - a ir )* m - W-'M , (A.12) 

[d^]a(l-y rQ )(— 1 )( 1 -)- 1 . (A.13) 

These equations are again just Eqs. ( |A.4|) and ( |A.5 ) with the replacement m — > m — 1. 



By the arguments of the previous section, we find it useful to modify the real-virtual 
counterterms similarly as done with the doubly-real singly unresolved terms. The only 
difference is the shift of m — » m — 1. The exact definitions are presented in Sect.[| and 
Sect-IH of the present paper. 

A. 3 Remarks on choosing do and d' Q 

As has been emphasized repeatedly, the integrals of the subtraction terms as defined in the 
present paper are m- independent for any do and d' . In this paper we set do = d' = 3 — 3e. 
Let us make a few comments on this choice. 

First of all, to avoid the introduction of spurious poles at oti r = 1 and/or y r Q = 1, we 
need to choose do and d' such that the overall powers of (1 — ai r ) and/or (1 — y r Q) that 
appear in any integral are non-negative. Although it is not manifest by considering only 
the singly-unresolved counterterms, it actually turns out that in the full NNLO scheme 
this implies 

do\ e=0 ,d Q \ €=Q >2. (A.14) 
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Secondly, we might want to avoid the appearance of negative powers of (1 — ati r ) and 
(1 — y r Q) in the subtraction terms themselves. This is because away from the limits both 
of these factors are between zero and one, thus if they are raised to a negative power, we 
are multiplying the subtraction terms by quantities that are greater than one, i.e. we are 
'over-subtracting' away from the limits. This leads us to choose 



*L 



=0 ,do| e=0 >m. 



(A.15) 



Since our primary interest is in 3-jet production in electron-positron annihilation, we set 
do\ e=0 = d' \ e=Q = 3. 

Finally, to fix the e-dependent part, we note that by choosing do = d' = m(l — e) the 
overall powers of (1 — ai r ) and (1 — y r g) in the doubly-real singly-unresolved subtraction 
terms (Eq. ( |3.2| ) and Eqs. ( |3.12| ) and fl3,13| )) are zero. Since we have already evaluated the 
integrated subtraction terms for this case (with ao = 1 and yo = 1) in Ref. [13], the choice 
of do = ^o = 3 — 3e is a natural one. 



B. Spin-averaged splitting kernels 



In this Appendix we recall the explicit expressions for the spin-averaged splitting kernels 
that enter Eqs. (O) and flOp. 



The azimuthally averaged Altarelli-Parisi splitting kernels read 
p£l( Zi ,z r ;e) = 2C A 



1 1 

\ £\ Z%Zr 

z% z T 



P q%r( Z ii Z r^) = T R 



- Zi z r 



p(%(z t ,z r ;e) = C F 
while their one-loop generalizations are 



1-e 



2 + (1 - e)z r 



(B.l) 
(B.2) 
(B.3) 



2C A 4 9 



l zeziZ r 

s. -, , ll Jifr = 99 , 

1 — e 



Pfcf T ( z ii z r'> e ) — r s'ren(^> z r\ e )-f/, ; / r ( z i-> z r '■> e ) + S , if fif r 

Cf r^ & (1 - ez r ) , if fif r = qg . 

(B.4) 
The renormalized rl^ n (z{, z r ; e) functions that appear above are expressed in terms of the 
corresponding unrenormalized ones as 



r S.,rcn\ z ii z r](-) — r s. \ z ii z r]^) „ r A \o 

le (47rj z cr 



/' 



2\ e 



cos 7re 



(B.5) 
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where the unrenormalized r^ r (zj, z r ; e) factors may be written in the following form 



33 ( ■ \ — A 

r S. \ z i: z r'i £) — o 

e z 



7re I Zi . , 

-r-, — r — +^ 2 i ? i(e,e, l + e,z r; 
sm(7re) \2 r> 



-^ e 2 Fi(-e,-e, 1 - e,z r ) 



(B.6) 



r a (.^i) ^r! £J 



-2 (C A - 2C F ) + - r 



vre / z, 



sin(7re) \z r 
ire { z 7 



+ 



l-2e 



sin(7re) \ z, 
A) — 3Cf 



+ C A - 2C F + 



+ z| 2 i ? i(e,e, 1 + e, z r 

-J +^ 2 -Pi(e,e, l + e,2j) 
Ca + 4T R (n f - n s ) 



3(3 - 2e) 



1 



2(C A - C F ) + Ca-^ ( - ) - CA«i 6 2-Fi(e, e, 1 + e, «r) 



-(C A - 2C F )z- e 2i ? i(-e, -e, 1 - e,*r) 



, (B.7) 



(B. 



The r^a non-singular factors are 

C A (1 - e) - 2T R (n f - n s 






(l-2e)(2-2e)(3-2e) ' 
For QCD, n s = of course. Finally /3q in Eqs. (|B.5| ) and 



„99 

' N.S. 



Ca - C F 
l-2e 



(B.9) 



is given in Eq. ( 4.15 ). 
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